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We investigate the Operator Product Expansion (OPE) on the lattice by directly measuring the 
product (J^Jv) (where J is the vector current) and comparing it with the expectation values of 
bilinear operators. This will determine the Wilson coefficients in the OPE from lattice data, and 
so give an alternative to the conventional methods of renormalising lattice structure function cal- 
culations. It could also give us access to higher twist quantities such as the longitudinal structure 
function Fl = F2 — 2xFi . We use overlap fermions because of their improved chiral properties, 
which reduces the number of possible operator mixing coefficients. 
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1. Introduction 

Our main theoretical tool for interpreting hadronic Deep Inelastic Scattering (DIS) is the Oper- 
ator Product Expansion (OPE). This relates the experimentally measurable electromagnetic tensor 

W^,{p,q) ^ {¥{p)\Ml)4i<l)\¥ip)) (1-1) 
to a sum over matrix elements of local operators. 

w,y{p,q) = {v{p)\Mq)JUq)\¥{p)) = Y,c';i,{q){Y{p)\ffnv{p)) ■ (i-2) 

m 

The local operators have interpretations in terms of the target hadron's internal structure. 

In each term of the OPE the scales separate. Dependence on the photon scale q is in the 
Wilson coefficient C'^y (q) while dependence on the quark momentum p is in the matrix element 
{Yip)\^'"\vip))- 

There has been a long history of lattice calculations of the hadronic matrix elements which 
occur in the OPE, but the Wilson coefficients are usually calculated perturbatively. In this work we 
investigate the possibility of also calculating the Wilson coefficients on the lattice by looking at the 
product of two electromagnetic currents. 

The Wilson coefficients are independent of the target. In our calculation we measure the 
current product (Compton amplitude) between quark states. We then plan to use the resulting 
Wilson coefficients together with lattice data on nucleon matrix elements to look at deep inelastic 
scattering. 

If we can measure the coefficients accurately enough we could learn something about higher 
twist effects, and non-leading amplitudes such as the longitudinal structure function Fl = F2 — 2xF\ . 
To calculate power corrections of this type we need to know both the matrix elements and the 
Wilson coefficients beyond perturbation theory 

In this work we report on an ongoing study of the lattice OPE using overlap fermions. These 
have the advantage of better chiral properties, which reduces the problems of operator mixing, and 
makes 0{a) improvement much simpler. Our earlier study using Wilson fermions was described 
in[§ 

2. Symmetry 

If p <^q SNQ will be able to truncate our set of operators according to their dimension. The 
operators we consider are the quark bilinears with up to three derivatives. These are the operators 
"VaFi//^, "V^FD^j \j7^D^^D^^ V^D^^D^^D^^ y where the matrix F can be any of the 16 matrices 
in the Clifford algebra. This means that there are a total of 16 x (1 + 4 + 4^ + 4^) = 1360 operators 
to consider. Do we need to find 1360 different C" values? 

To reduce the number of independent coefficients we want to choose a q vector with as much 
lattice symmetry as possible, so we have taken qoc (1,1,1,1). The data presented here are for the 
choice aq= (f ,f ,f ,f)- 

In the expansion for W44 we know that rotations and reflections that mix the space direc- 
tion x,y,z are symmetries. So it is obvious that the Wilson coefficients of l/AyiDii/A, YY2D2W 
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xffYiD^Y should all be the same. Exploiting symmetry this way reduces the original 1360 coeffi- 
cients down to only 67, a much more manageable problem. We can also use symmetries to relate 
different components of W^v For example, the Wilson coefficient of 1/774D4I/A in W44 is the same 
as the C for 1/77303 i/a in W33. 



3. Lattice Details 

We are carrying out calculations with overlap valence fermions, due to their superior chiral 
properties. However because of the high cost of a full dynamical overlap calculation, we have to 
use gauge configurations calculated with Nf = 2 dynamical clover fermions. 

The overlap fermions are calculated with p = 1.5 and bare mass am = 0.024. The results 
discussed in this work used a 16^ x 32 lattice, simulated with Nf = 2 clover fermions, at j8 = 5.29 
and K = 0.1350, which corresponds to a lattice spacing of a = 0.075 fm. 

The Green's functions have been 0{a) improved, which is easy to do with overlap fermions. 
The current /^(jc) was represented by the local current 1/7(;c)7^i/a(x). Since we are measuring op- 
erators between quark states we have to fix the gauge. We used the lattice Landau gauge. We use 
a momentum source [||, 0] for all Greens functions, which leads to a great reduction in statistical 
noise. 

We used a momentum transfer aq = (fifjfjf) corresponding to |<7| = S.3GeV. For this 
q value we measured the two-point and three-point functions for 28 different p vectors, with a 
large spread in directions. Here we give results for {JaJa), we are also collecting data on other 
components of the current-current tensor. 

We only consider the flavour non-singlet case, so we do not include any purely gluonic opera- 
tors in our calculation. 



4. Strategy 



We calculate the Compton scattering amplitude for a quark with a large number of p values. 
The result is a Dirac matrix for each p vector, i.e. 16 complex numbers. So our data on the Compton 
amplitude consist of 16Mp complex numbers. We also calculate the operator Greens functions for 
each of our operators, for all of these p values. This gives (16Mp) x No numbers as our data on the 
operators. 

The information we want to extract from all this data are the No Wilson coefficients which 
best reproduce the Compton scattering amplitudes. 

This is essentially a linear algebra problem, and is best written as a matrix equation 



'^iiv{Pnq)=Y.{w{Pi)W"'W{Pi))C;,{q) 
m 



(16Mp) 




(4.1) 



(4.2) 



{l6MpXNo)-{No 
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5. Singular Value Decomposition (SVD) 

There are two difficulties with this system of equations. Firstly if 16Mp > No the system is 
overdetermined, so no exact solution will be possible. Secondly, some of the operators might be 
linearly dependent, in which case the system is also ill-conditioned, and there will not even be a 
unique best approximation to the solution. 

Nevertheless by using Singular Value Decomposition (SVD) [Q], the standard technique for 
problems of this type, we can find values of the coefficients C, which give a very good, and very 
stable, approximate solution. 

We can always factorise the operator matrix O as 

O = UcoV^ (5.1) 

N 

Oij = £ UikCOkVjk (5.2) 

where U and V are column-orthonormal, U is (16Mp) x No, V is No x No and ft) is a diagonal 
No X No matrix, with positive real eigenvalues cot arranged in descending order. The (Ok are the 
analogues of eigenvalues for a rectangular matrix, while the matrices U and V contain the eigen- 
vectors for the system. If some of the (Ot are very small they can safely be dropped from the sum 
without significantly changing O. 

We define O^"^ as the approximation to O that we get by keeping the n largest (O values, and 
dropping the N — n smallest values. This gives 

Otj = I UikCOtVjt = f/(")ft)(")(v(«))^ (5.3) 

where now is (16Mp) x n, ftj(") isnxn and V^"^ is No x n. 0^") is still (16Mp) x No- It can 
be shown that the least-squares solution to 0^"^C = WisC = V^"\(o'-"^y\U^"^yW . 

What singular value decomposition has done is to find the N — n linear combinations of our 
original operators which have the least influence on W, and discarded those combinations. Dis- 
carding these operators makes the problem much more stable. However if we discard too many 
operators, O^"' will no longer be a good approximation to O, and so the solution of the approxi- 
mate system will not give a good approximate solution to the true problem. 

To solve a system by SVD we vary n, the number of COt values retained, looking for a region 
where the residue R = {W — OC)^ is small and the system is stable. Fig.|I] shows how the residue 
declines as the number of singular values is increased. A plateau is reached beyond n Rr: 40. Beyond 
this point adding more operators does not decrease the residue significantly. 

To judge the stability of the fit, we look at the value of one of the Wilson coefficients, and 
see how it depends on n. In Fig.^ we show the result for the choice C2 in our list of operators, the 
operator \j7{YiDi + 72D2 + 7303)1/^. At first the coefficient changes dramatically as operators are 
added, but by the time n has reached 40, there are only minor changes. If n is made too large, there 
is a risk that we will start "fitting to the noise", and the value of the coefficient will become noisy. 
There are indeed some fluctuations beyond n = 50, but they are not unduly large. 

Another way of judging the quality of the fit is to exclude one momentum value from the fit, 
and see how well it is predicted by the data at all the other momentum values. The result is shown 
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Residue vs singular values 
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Figure 1: The residue of the fit as a function of n, the number of singular values retained. 

Wilson coefficient vs singular values 

0.035 I , , , , , , 1 




n 

Figure 2: The value of C2, the Wilson coefficient for '^f{y\D\ + Y2D2 + 73/33)'/, as a function of n, the 
number of singular values retained. 

in Fig.^ Again we see that we need at least 40 singular values to produce a fit with good predictive 
power. If, at large n we started fitting to the errors in the data, we would expect to see the predictive 
power becoming worse. There doesn't seem to be any sign of this happening. 

6. Results 

In Fig.^ we show the results of our fit. Chiral symmetry shows up well. Operator number 1 , 
the operator VAlyA, and operators 7 to 16, which are two-derivative operators proportional to the unit 
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Figure 3: The accuracy with which the fit predicts the value of W^v as a function of n. 

matrix or to a matrices, are ruled out by chiral symmetry, so they would have Wilson coefficients 
of in the chiral limit. Here, using overlap fermions, we find that these coefficients are indeed 
small. This contrasts with earlier work with Wilson fermions, where the operator 1/71 

was prominent. 




Figure 4: The Wilson coefficients for determined on the lattice according to our procedure. They are 
compared with the lattice tree-level values, shown by the blue line. 

We have calculated the tree-level Wilson coefficients for overlap fermions, following the cal- 
culation set out in [^. The tree-level results are shown by the blue line in Fig.^ We see that the 
pattern of Wilson coefficients is very similar. 
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These results are at a single value of We intend to gather data at other values to investigate 
the size of lattice artefacts, which should be 0{cP') for overlap fermions. 

7. Conclusions 

Using Singular Value Decomposition we have been able to reconstruct Wilson coefficients 
from lattice data on the electromagnetic tensor. Statistical errors are small, due to our use of 
momentum sources for the inversions. 

The results we have look reasonable, they follow a pattern similar to that seen at tree-level, 
and they show the effects expected from chiral symmetry. So far we only have data from a rather 
large value of |a^|, we plan to look at more q values to check for 0{a^) lattice errors. 
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